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We investigate the properties of dense matter and neutron stars. In particular we discuss 
model calculations based on the parity doublet picture of hadronic chiral symmetry. In 
this ansatz the onset of chiral symmetry restoration is refiected by the degeneracy of 
baryons and their parity partners. In this approach we also incorporate quarks as degrees 
of freedom to be able to study hybrid stars. 

1. Introduction 

One of the main and still unresolved problems in nuclear physics is the understand- 
ing of the properties and phase structure of hot and dense matter. Here heavy-ion 
physics is an important experimental approach to extract information on very hot 
and, depending on the beam energy also dense matter. As a complementary ap- 
proach the investigation of neutron stars can yield constraints on the properties 
of extremely dense and rather cold matter, where the temperature T is usually 
less than 1 MeV except for the short period of the proto-neutron star, where T 
might reach 30 to 50 MeV. Therefore both fields together can help to pin down, 
or at least constrain the equation of state and other properties of strongly inter- 
acting matter far beyond the nuclear matter groundstate. One important aspect 
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of QCD is its inherent chiral symmetry that characterizes the transition from low- 
density /temperature matter to matter at high values of T and p crossing from the 
symmetry broken to the restored phase. There are different ways to construct a chi- 
rally symmetric hadronic model. One approach, given by the linear <j model and its 
many variations and extensions. Here the masses of the baryons are largely gener- 
ated by the coupling of baryons to the scalar field. Therefore, in the chirally restored 
phase the baryons become very light. Alternatively, in the so-called parity-doublet 
approach the baryonic degrees of freedom are augmented by the opposite-parity 
counterparts. ^ One possible state of such a particle is the N(1535) as parity part- 
mer of the nucleon. In this formulation it is possible to define a mass term for the 
baryons without violating chiral symmetry. The chirally restored phase is charac- 
terized by the degeneracy of both parity partners. It has been sh own tha t it is 
possible to describe saturating nuclear matter in such an approach .I^EllEI in the 
following we discuss results based on the latter approach, studying hadronic SU(2) 
and SU(3) models as well as an extended description including quark degrees of 
freedom. We present results for the phase structure of the models as well as the 
properties of compact stars. 

2. The Hadronic SU(2) Parity Model 

In the parity doublet mo del positive and negative parity states of the baryons are 
arranged in doublets The two components of the fields defining the respective 
parity partners, and (f- transform in opposite way regarding chiral transforma- 
tions: 

if'_^ = Lif-R ,(p-L=R^-L' (1) 

Here, L and R denote the rotations in the left- and right handed subspaces. This 
allows for a chirally invariant mass term in the Lagrangian of the general form: 

^o(^-75^+ - ^+75^-) = 

mo{(^-L^^R - (p-R^+L - (P+L^-R + (p+R^-h), (2) 

where mo represents a mass parameter. The coupling terms of the scalar field to 
the nucleonic states read 

Cbs = [cr + Z75r + hCp- {a - ij^r • tt) . (3) 

As both states in the doublet mix one has to diagonalize the states as shown in^. 
Fixing the coupling constants to reproduce the vacuum masses of the nucleons the 
resulting effective masses of the diagonalized states read ^ : 



(Mjv+ + M^_)2 ■ 
i "^0 
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where ctq denotes the vacuum value of the scalar field. As one can see from this 
expression both states are degenerate in the limit of a vanishing scalar field. 

Doing an exhaustive parameter scan and constraining the saturation density po 
to be about 0.15//m^ and the binding energy per nucleon B/A to be around 16 
MeV the resulting compressibilities k of nuclear matter are shown in Figure 1 as 
function of the mass parameter mo-^ One can observe a strong dependence on this 
parameter. In order to have a reasonable value of n below 300 MeV a rather high 
value of mo around 840 MeV has to be chosen. 




Fig. 1. Left Panel: Compressibility tx as function of mass parameter mo given the mass of the 
nucleonic parity partmer Mjv- = 1535MeV. Right Panel: Corresponding maximum neutron star 
masses in units of the solar mass. 



Using the precious parameter scan with the same constraints on po and B/A 
and additionally demanding that k has a value below 400 MeV the Tolman- 
Oppenheimer-Volkoff equation is solved for all acceptable parameter sets, subse- 
quently determining the maximum star masses. The result of this calculation can 
be seen in Figure 1 (right panel). 

As the identification of the nucleonic parity partner is unclear one can study 
the model for different values of Mat-. Below a value of about 1380 MeV the parity 
partner states get populated in the core of the neutron star. 

Setting the mass value to 1200 MeV (Fig. 2) a substantial fraction of the baryons 
inside the star are given by the parity partners that start to appear at about 2po in 
the case of n_ and 3po for p_.^ One can observe that for large baryon densities the 
respective densities of the nucleons and the parity partners become equal, which is 
a consequence of the degeneracy of the two states in the chirally restored phase. 

This restoration of chiral symmetry can be seen in the behavior of the scalar 
field as it is shown in the right panel of Fig. 2. Here also the difference of the 
drop of the scalar field value in the case of isospin symmetric and star matter is 
shown. As it is generally the case the transition becomes weaker and shifts to lower 
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Fig. 2. Left Panel: Normalized particle densities as function of baryon density for star matter with 
a mass of the nucleonic parity partner of Mjv- = 1200MeV (parameter set PI in Ref. ^). Right 
Panel: Scalar field a as function of density 



densities for star matter. This is due to the breaking of the isospin symmetry in the 
stehar environment. Thus, the parity partner of the neutron is populated earher, 
but as this onset of population affects only the n_ compared to the simultaneous 
appearance of the n_ and p_ in the case of isospin symmetric matter the effect on 
the scalar field is reduced. 



3. The SU(3) Parity Model 

Extending the model to flavor SU(3) in a nonlinear realization of chiral symmetry 
we follow the discussion in In general the formalism is analogous to the SU(2) 
case (see Ref. ^ for a detailed discussion). The baryonic Lagrangian reads 

i i i 

with couplings to the vector fields and (j). The effective masses of the baryon 
octet and the corresponding parity partners are given by the expression 



m 



i± 



including a SU(3) symmetry breaking term proportional to the number of strange 
quarks of the baryon. rris has been fixed to a value of 150 MeV. The index i 
denotes the baryon and ± are the positive and negative parity partners as in the 
SU(2) case. Note that in the SU(3) formulation there are two scalar fields a and 
which are effectively the counterparts to the QCD scalar non-strange and strange 
quark condensates. The scalar meson interaction driving the spontaneous breaking 
of chiral symmetry can be written in terms of SU(3) invariants Ii = Tr(I]) , I2 = 
Tr{T?) , h = det{^) , h = Tr{T.'^), where S denotes the multiplet of scalar 
mesons?^ 

V = Vo + \k0I2 - hii - k2h + kMh) , (7) 
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where Vq is fixed by demanding a vanishing potential in the vacuum. The exphcit 
symmetry breaking term that generates the correct pion and kaon masses with their 
corresponding decay constants can be written as ^ 

LsB = mlUa + (y2mlfk - 75^^ A) C- 

(8) 

Fitting the couphng constants to a reasonable nuclear matter ground state as 
in the previous section we solve the equations of motion for the various fields. In 
a further extension we include quark degrees of freedom a nd th e Polyakov loop 
field ^ as order parameter for the deconfinement transition EUHl Here we follow 
the formulation as outlined in Ref. The quarks couple to the condensates and 
attain an effective mass: 

ml=gq^a^moq, m*=^5cC + ^0s, (9) 

with couplings gq^ = QsC =4.0. moq and mos describe the current quark masses 
as well as contributions from the slowly varying gluon condensate The potential 
of the Polyakov loop effectively represents the contribution of the gluons to the 
thermodynamical quantities of the system like pressure and energy density. Here 
we adopt the ansatz proposed in^. 

U = -ia(T)^^* + b{T)ln[l - 6<I><I>* + 4(<I>^$*^) - 3(<I><I>*)2] (10) 

with a(T) = aoT"^ + aiToT^ + asTo^T^, b{T) = h^T^T, where Tq sets the scale for 
the first-order deconfinement phase transition in the quenched limit o f only gluons. 
Finally, we introduce an excluded volume correction for the hadrons which 
leads to the suppression of hadronic degrees of freedom relative to quarks at higher 
densites as described in^. 

Using the parameter set for a mass of the negative parity state of the nucleon 
Mat- = 1535MeV and assuming, for simplicity, the same mass splitting of all 
baryon doublets as the nucleonic one, as the possible attribution of observed states 
is less clear for the hyperons, we study the phase transition behavior of the model 
as function of chemical potential and temperature. The result is shown in Fig. 
4. The two plots are calculated for different values of the temperature parameter 
in the Polyakov loop potential Tq = 270 and 220 MeV, respectively. The lower 
band indicates the chiral transition region during which the value of cr drops from 
0.8 to 0.2 of its vacuum value. The upper band corresponds to the deconfinement 
transition connected to the value ^ of the Polyakov loop changing from 0.2 to 0.8. 
The two first-order transition lines for the liquid-gas and chiral transition are shown 
as solid lines. Both stop at a critical end-point with values (in units of MeV) of 
(T,/i) = (17,903) and (58,1200), respectively. 

In order to study the particle species that are populated the densities of the 
various degrees of freedom as function of chemical potential are shown in Fig.lH The 
results in the left panel are obtained for zero temperature and isospin-symmetric 
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Fig. 3. Phase transition diagram as function of temperature T and quark chemical potential 
= M-B /3. The two panels are the results for different parameter values for To of 270 (left) and 
220 MeV (right). The lower band shows the range where the scalar field has a value between 20 and 
80 percent of its vacuum value. Analogously, the upper band indicates the range of temperatures 
and chemical potentials for which the Polyakov loop field lies between 0.2 and 0.8. The solid lines 
show the first-order liquid-gas and chiral phase transitions. In case of the lower value of To the 
baryochemical potential of the critical end point moves from 1200 to 1150 MeV. 

matter. One can see the onset of the occurrence of quarks and the nucleonic parity 
partners at shghtly higher chemical potential, such that there actually is a two-step 
first-order transition. In contrast to the hadronic model the parity partner appears 
earlier driven by the appearance of the quarks and the connected drop of the scalar 
condensate. A similar, but more involved behavior can be seen in the right panel 
for star matter. First the down quark appears in a first-order transition, then, in a 
second transition the positively charged up quark and the nucleonic parity partners 
are populated. The hyperonic states and their parity partners are suppressed by 
the occurrence of the quarks and appear only at much higher chemical potential. 

Using the thus obtained equation of state we solve the TOV equations varying 
the central pressure of the star. The resulting star sequences are given in Figure 
[5l The result for the purely hadronic model with a maximum mass of about 1.6 
solar masses and a radius of 12 km. Including quarks shifts the maximum mass to 
1.96 solar masses with a reduced radius of 10.7 km, which is in agreement with the 
rather accurate measurement of the 1.97 ± 0.4 solar mass pulsar PSR J1614-2230. 

The kink in this sequence indicates the onset of a quark core in the star, i.e. 
the transition from a neutron to a hybrid star. The fact that the masses for both 
cases below that value are not identical originates from the fact that only in the 
case of the full calculation with quarks the excluded volume corrections are taken 
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Fig. 4. Normalized baryon number densities for the different particles as function of chemical 
potential for isospin symmetric matter (left panel) and beta-equlibrated stellar matter (right 
panel) at T = 0. 

into account. 
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Fig. 5. Mass-radius diagram for compact stars in the SU(3) parity model including hadrons (h), 
and quarks and hadrons (h+q). 

Summarizing, we discussed SU(2) and SU(3) parity doublet models for hadronic 
matter with and without quark degrees of freedom. We could reproduce nuclear 
matter saturation properties reasonably well. Depending on the mass of the parity 
partners the model has a first-order chiral phase transition at high densities for 
larger masses and a crossover for smaller masses, respectively. Including quark de- 
grees of freedom triggers a first-order transition for temperatures below ~ 60 MeV. 
The calculation of static hybrid stars yields maximum star masses of 2 solar masses 
in agreement with observation. 
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